Low-temperature specific heat and thermal conductivity of glycerol 
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We have measured the thermal conductivity of glassy glycerol between 1.5 K and 100 K, as well as 
the specific heat of both glassy and crystalline phases of glycerol between 0.5 K and 25 K. We discuss 
both low-temperature properties of this typical molecular glass in terms of the soft-potential model. 
Our finding of an excellent agreement between its predictions and experimental data for these two 
independent measurements constitutes a robust proof of the capabilities of the soft-potential model 
to account for the low-temperature properties of glasses in a wide temperature range. 

PACS numbers: PACS: 65.40.+g, 66.70.-|-f, 63.50.-l-x 
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It is well establishedH'EI that glasses or amorphous solids 
exhibit thermal properties very different from those of 
crystalline solids and, even more strikingly, very similar 
among themselves irrespective of the type of material, 
chemical bonding, etc. At temperatures T < 1 K, the 
specific heat Cp of non-metallic glasses is significantly 
larger and the thermal conductivity k orders of magni- 
tude lower than those found in their crystalline coun- 
terparts. Cp depends approximately linearly (Cp oc T) 
and K almost quadratically (k cx T^) on temperature, in 
contrast to the cubic dependences observed in crystals 
for both properties, well understood in terms of Debye's 
theory. At T > 1 K, Cp still deviates strongly from the 
expected CDobyc oc dependence, exhibiting a hump 
in Cp/T^. In the same temperature range the thermal 
conductivity exhibits a universal plateau. 

In lfl72, Philfipfi and Anderson, Halperin and 
Varmaa introduced independently the well-known tun- 
neling model (TM), whose fundamental postulate was 
the ubiquitous existence of small groups of atoms in 
amorphous solids which can tunnel between two con- 
figurations of very similar energy. This simple model 
of two-level systems or tunneling states successfully ex- 
plained the low-temperature properties of amorphous 
solidsO, though only for T < 1 K. On the contrary, the 
also rich universal behavior of glasses above 1 K (the 
hump in Cp/T^ and the plateau in the thermal con- 
ductivity, or the remarkable feature in the vibrational 
density of states g{v)/v^ at low frequencies known as 
the boson peak) still remains a matter of debate. One 
of the best accepted approaches to understand all the 
general behavior of glasses in the whole range of low- 
energy excitations is the phenomenological soft-potential 
model (SPM), which can be regarded as an extension 
of the TM. The SPMeIlI postulates the coexistence of 
extended lattice vibrations (sound waves) with quasilo- 
calized low-frequency (so/t) modes. In this model, the 
potential of these soft modes has a uniform stabilizing 
fourth-order term W. In addition, each mode has its 
individual first-order asymmetry Di and second-order 
restoring force terms D2, which can be either positive 
or negative. Similarly to the TM, a random distribution 



of potentials is a^med^(I)i, D2) = A-|The SPM has 
been developedaQu'ElEEfEil and reviewedtSllj in earlier pa- 
pers where the interested reader is referred to. 

Glycecpl [€3115(011)3] is probably the most widely 
studiedtil glass-forming liquid. Its high viscosity at a 
melting point around room temperature (T^ = 291 K) 
provides experimentalists with a very convenient temper- 
ature range where the supercooled liquid can be stud- 
ied. Below the glass transition at Tg ~ 185 K, the 
frozen-in liquid becomes a glass with a relatively weak, 
hydrogen-bonded network structure. Despite its good 
glass-forming ability, glycerol Gap also be obtained in an 
orthorhombic crystalline statell3, with four C3H5(OH)3 
molecules per unit cell, bHulding up a structure of infinite 
hydrogen-bonded chaingla. Several measurements of the 
specific hea,t || 0f | t his well-known glass have been indeed 
reportedllIrll3ll3l23, though not reaching temperatures be- 
low 1.5 K in an y c as e. A broad maximum in Cp/T^ was 
clearly observecOE^I around 8.5 K, but the expected exis- 
tence of tunneling states could not be determined, since 
it requires temperatures, typically below 1 K. Further- 
more, Calemczuk et aLtHI also measured the specific heat 
of the crystalline state of glycerol, but only down to 5 K, 
hence not reaching temperatures low enough as to assess 
its Debye coefficient. On the other hand, there are no 
published thermal-conductivity data at low temperatures 
neither for glycerol nor for any other similar molecular 
glass. This is very probably due to the experimental dif- 
ficulties of adapting the thermal-conductivity technique 
to a sample which is liquid at ambient temperature and 
must be thermally controlled in situ to freeze it into the 
glass state, maintaining at the same time an appropriate 
geometry and a moderate heat flow for the thermal con- 
ductivity to be correctly measured. To our knowledge, 
only a few orientationally-disordered ("glassy") cm^als 
from other molecular liquids have been measuredcilEHI. 

In this work, we report thermal-conductivity data of 
glassy glycerol measured from 1.5 K to around 100 K, 
as well as specific-heat measurements of both glassy and 
crystalline glycerol between 0.5 K and 25 K. In addition, 
we make concurrent use of both thermal properties mea- 
sured for the glassy state in order to test the validity of 



2 



the soft-potential model in a typical molecular glass. 

In Fig. |l|(a), we show our thermal-conductivity data 
of glassy glycerol obtained by cooling in situ below Tg 
at a rate around -1 K/min liquid glycerol (Merck, anhy- 
drous, used without futher purification), placed inside a 
very thin-walled (0.2 mm) nylon tube. Standard steady- 
state techniques were employed. The thermal conductiv- 
ity of the empty nylon tube was independently measured 
and subtracted. As can be seen, glycerol exhibits the 
thermal-conductivity behavior typical of glasses, with a 
plateau around 10 K quantitatively very siniilar to that 
of strong, network glasses as AS2S3 and Ge02H, i.e., a rel- 
atively high thermal conductivity among glasses. There- 
fore, glycerol could be a useful heat exchange medium at 
low temperatures when used as a glassy matrix. 

Let us first analyze the thermal conductivity of glyc- 
erol within the soft-potential model. According to the 
SPMQ, the inverse mean-free path of the phonons car- 
rying out the heat is simply the sum of three contribu- 
tions: the resonant scattering of the sound waves by ei- 
ther tunneling states or by soft vibrational modes, and 
the scattering by classical relaxational processes in the 
same asymmetric double-well potentials responsible of 
the tunneling states. From these premises, and insert- 
ing the corresponding SPM expressions for these three 
inverse meaitfc^e paths, the thermal conductivity k was 
found to hM 
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F{z), 
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where 



F{z) 



dx- — 
(1 



-^)2 1.1 tanh(2;/2) + 0.7z3/4 ^_ ^3^3^ 

_ (2) 
C is the usual (averaged over longitudinal and trans- 
verse acoustic modes) dimensionless constant of the TM, 
directly related to the universal plateau in the internal 
friction = (7r/2)C, and z = k^T/W, where the en- 
ergy W is the aforementioned parameter of the SPM, 
which marks the crossover from the tunneling-states re- 
gion at the lowest temperatures to the soft-mades region 
above it. Indeed, W can be determinedllinij from the 
position of the maximum Tmax in a n/T versus T plot: 



W ~ l.QkBT„ 



(3) 



Tmax therefore separates the low-temperature range 
where resonant scattering by tunneling states domi- 
nates (k oc T^) from higher temperatures {k « const) 
where the soft modes are the main scatterers of acoustic 
phonons. 

Fig. |l| shows that SPM eqs. provide a good fit 

to thermal-conductivity data of glycerol from the lowest 
temperatures up to the end of the plateau around 3|0-|K, 
The two parameters, W and C are easily obtainedM't^ 
from the k/T plot [Fig. |l](b)]: the maximum position 
gives IF = 4.3 K and its height C = 1.9x10"''. The latter 
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FIG. 1: Low-temperature thermal conductivity of glassy glyc- 
erol, (a) : Thermal conductivity as a function of temperature 
in a log- log plot. Solid line is a fit to the soft-potential model 
(SPM). (b): Thermal conductivity divided by temperature 
below 10 K. The maximum marks the crossover from tunnel- 
ing states to quasilocalized vibrations (soft modes) as domi- 
nant scattering centers for the "phonons" and allows a direct 
determination of the SPM parameter W . 



constitutes indeed a prediction for the low-temperature 
internal friction of glycerol, never measured to our knowl- 
edge: Jisjplateau value should be around Q"^ « 3 x lO"**, 
up toEltJ 1.2W/ fee ~ 5 K, where the rise due to thermal 
relaxation should occur. 

Our specific-heat measurements of glassy and crys- 
talline states of glycerol are presented in Fig. || (in a 
log-log plot showing the whole low-temperature range) 
and in Fig. ^ (in a Cp/T vs plot at the lowest tem- 
peratures to address the tunneling-states range for the 
glass and the Debye limit for the crystal). The heat ca- 
pacity was measured in a ^He-cryostat, employing a low- 
temperature quasi-adabiatic calorimetiiijipll, similar to 
one previously used in a ^He-cryostatcSEj. The liquid 
glycerol is placed in a vacuum-tight, thin- walled copper 
can, with a fine mesh of copper fitted inside to facilitate 
thermal equilibrium. A thin gold wire was used as heat 
switch to cool the experimental cell. The subtracted ad- 
denda contibution to the total measured heat capacity 
at 4.2 K (1 K) was about 15 % (22 %) for the glass and 
about 30 % (50 %) for the crystal. The glass state was 
obtained by simply cooling the liquid from room temper- 
ature down to liquid- helium temperatures. After having 
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measured its heat capacity at low temperature, glassy 
glycerol was slowly heated above Tg (calorimetrically ob- 
served to take place at 185 K) until it crystallized around 
250-260 K. Once this first-order, exothermic transition 
was completed, the sample was cooled again and the heat 
capacity of the crystal measured. 
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FIG. 2: Low-temperature specific heat C-pjT'^ of glass (Q) 
and crystal (□) phases of glycerol for the whole measured 
temperature range in a log-log scale. Solid triangles are pub- 
lished data from Leadbetter and WycherleytJ. Dotted lines 
show the correspondent Debye contributions to the specific 
heat either measured calorimetrically (crystal) or estimated 
from sound velocities (glass). Solid line shows the curve cal- 
culated with the soft-potential model (SPM), taking the pa- 
rameter W from the thermal conductivity data. Dashed line 
shows the same SPM calculation, but with a gaussian cutoff 
of the asymmetry for soft modes at high energies (see text). 
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FIG. 3: Low-temperature specific heat of glass and crystal 
phases of glycerol in a Cp/T vs plot. Symbols are as in 
Fig. ^. Solid lines are least-squares linear fits. 

In Fig. 1^, published data for the glass-between 1.9 and 
25 K from Leadbetter and Wycherleyll^ can be seen to 
show a very good agreement with our data. There is also 
a good agreement with data for the crystal above 5 Kl3. 
The low-temperature cubic Debye contributions are also 
indicated in that figure: for the crystal, it has been ob- 
tained from a least-squares linear fit below 1.5 K in Fig.|^, 
giving a Debye temperature 8d — 367 K; for the glass, 
it can be estimated from zero-temperature extrapolations 



of elastic dataSBEto be Od = 317 K. In contrast to the 
crystal, where Cp holds the expected Debye behavior for 
T < 0d/5O, the glass exhibits the typical excess over the 
Debye contribution in the whole low-temperature range. 
Within the SPM, these additional excitations are under- 
stood as soft quasiharmonic vibrations in single-well po- 
tentials responsible for the "boson peak" in Cp/T^, to- 
gether with tunneling states arising from related double- 
well potentials responsible for the quasilinear contribu- 
tion. The latter dominates below thejaaiaimum Tmin in 
Cp/T^ given by « 1.8 - 2 fcBTminBEl, and is more 
clearly observed in Fig. ^ where a linear coefficient 7 — 
0.13 mJ/mol-K^ can be obtained (i.e., 7 = 1.4^J/g-K^ 
that is very similar to those values found in typical net- 
work glassesEloEl). 

Finally, we wish to compare the specific-heat data of 
glassy glycerol with the behavior predicted by the SPM, 
once we have previously obtained its basic parameter 
W — 4.3 K from thermal-conductivity data. As said 
above, the SPM postulates the coexistence of extended 
sound waves with quasilocalized modes, either tunneling 
states in double-well potentials or quasiharmonic vibra- 
tions in single-well potentials, with a gradual crossover 
between them characterized by the energy W . There- 
fore, we can write the specific heat for a glass as 

Cp — CDcbyc + CtLS + C'sm, (4) 

where Cgm is the contribution of soft modes (see below), 
and the contribution of the tunneling states Ctls is de- 
termined, as in the TM, by m^ans of the well-known ex- 
pression for two-level systemsB, using the correspondent 
densitv-of tunneling states in terms of the SPM (eq. (4.7) 
in RefO or eq. (9.30) in Ref.B). 

In the standard SPM, the density of quasiharmonic 
soft vibrations increases continuously with frequency as 

ff.n.(H-^§(^) , (5) 

where is the distribution constant of soft potentials. 
In the harmonic approximation, this leads Xq a specific- 
heat contribution from soft modes given byli3 

a„ = |!p...e(^sI)^ „) 

The solid line in Fig. ^ shows the result of eq. 
(^), where W = 4.3 K was taken from the thermal- 
conductivity measurement and Pg — 1.6 x 10"'^^ mol"""^ 
was simply determined to scale with experimental data. 

Obviously, the real distribution of soft modes cannoJ. 
grow with frequency gs-m{v) oc v'^ unlimitedly. Gil et als 
proposed a gaussian distribution in the asymmetry of the 
soft potentials (hence multiplying eq. (||) by an integral 
factor, see eq. (9.40) in Ref.El) based in a thermal strain 
ansatz, which without any further fitting parameter al- 
lowed them to account for the specific heat, thermal con- 
ductivity and vibrational density of states g{v)/v^ in the 
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whole relevant range, inplpding the "boson peak" . Alter- 
natively, Gurevich et aZ.E3 argued that the simple picture 
of independent quasilocalized harmonic vibrations coex- 
isting with sound waves should fail at the loffe-Regel 
limit, not far above the boson peak. The interaction 
between soft modes would lead to a reconstruction of 
the vibrational density of states at higher frequencies, 
where delocalized soft vibrations with gsm('^) oc v should 
dominate, therefore also explaining the boson-peak fea- 
ture. Furthermore, hybridization of acoustic phonons 
with quasilocalized modes has been proposed to set in 
around the boson peakcS. For the sake of completeness, 
we also show in Fig. ^ by a dashed line the result of the 
total SPM predictioit-fior the specific heat if the above- 
mentioned correctionBE^ of eq. (^) is used. As can be 
seen, the position of the maximum in Cp/T'^ (which only 
depends on W and the glass transition temperature Tg, 
hence being determined without any free parameter) is 
very well predicted, though its height is not so well ac- 
counted for. However, these quantitative agreements or 
disagreements of the SPM around or above the boson 
peak are perhaps not very relevant, since the low-energy 
limits of both independent quasilocalized modes and De- 
bye acoustic phonons (note that even in the crystal, Cp 
starts to deviate from the cubic limit above 5 K) should 
begin to fail there, for the reasons mentioned above. Nev- 
ertheless, one may conclude that the maximum in Cp/T^ 
of glasses (i.e., the boson peak) is just the fingerprint 
of the end for the low-energy distribution of independent 



soft modes, which govern low-temperature properties and 
low- frequency dynamics of glasses. 

In any case, it is noteworthy that with only a con- 
stant factor Pg, used as free parameter to fit the height 
of the Cp{T) curve iW was independently determined 
from thermal-conductivity data), the SPM is able to ac- 
count consistently for the specific heat of this glass in the 
low-tcmperature range, from the tunneling- states region 
below I K up to the broad peak in Cp/T^ , including the 
crossover region around the minimum in Cp/T^ . 

In summary, we have concurrently measured two low- 
temperature thermal properties of a paradigmatic molec- 
ular glass, glycerol, as well as the specific heat of the 
crystalline phase between 0.5 K and 25 K, hence being 
able to determine its Debye temperature. At lower tem- 
peratures, both properties for the glass exhibit a typical 
behavior indicative of the existence of tunneling states. 
Moreover, we have used these data as a new test of the 
SPM, which has been shown to successfully explain the 
specific heat and the thermal conductivity in a wide tem- 
perature range, also for a molecular glass such as glycerol. 
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